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I. INTRODUCTION 

Scalar fields occupy a central place in modern cosniol- 
ogy.On one hand, these fields (or field) play the role of 
infiatons, i.e., fields responsible for the inflation in the 
early universe 0-0 ■ On the other hand, massless scalar 
fields are currently considered a quintessence, a physi- 
cal substratum responsible for the observed accelerated 
expansion of the universe > ■ 

According to the most advanced cosmological model 
today, the so-called chaotic infiation theory it is as- 
sumed that the early universe was dominated by a scalar 
field with minimal coupling (in the simplest case) and 
that the energy density was concentrated in the self- 
action potential. Under these conditions, the initial equa- 
tions can be considerably simplified. After the simplifi- 
cation, the obtained system is easily integrated, showing 
the existence of an inflationary phase (more accurately, a 
quasi-de Sitter phase during which the Hubble parame- 
ter can be considered constant). It is then assumed that 
during the inflation process, the kinetic term increases 
until the slow rolling approximation is no longer appli- 
cable. It then follows that inflation ends spontaneously. 
The next phase is the oscillation phase, which is neces- 
sary to flU the universe, "emptied by inflation," with ele- 
mentary particles. According to the discussed paradigm, 
all these particles were produced from the vacuum by a 
rapidly oscillating scalar field whose oscillation amplitude 
decreased with time. A secondary reheating occurred, af- 
ter which the universe not only was homogenous and flat 
but also was filled with hot matter. In other words, at 
this point, we can use Friedmann's cosmological equa- 
tions for matter with a state equation characteristic of 
an electromagnetic field. The subsequent evolution of 
the universe follows the classical scenario: the universe 
cools as it expands, radiation decouples from matter at 
some point, and the universe then expands according to 
Friedmann's "two-thirds" law, and so on. 
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The discovery of the current accelerated expansion of 
the universe [4] was unexpected because this fact did not 
follow at all from the existing inflation paradigm. The ex- 
planation demanded additional hypotheses and assump- 
tions about a quintessence (associated with a "dark en- 
ergy," which is again "overabundant" in the universe) 
that is the cause of the universe's accelerated expansion. 
The need for additional assumptions of course contra- 
dicts the established scenario. Perhaps, the best that 
can be done to preserve the current paradigm is to prove 
that both the early inflation and the current acceleration 
are caused by the same phenomenon, i.e., the same scalar 
fleld. If this proves impossible, then the inflation scenario 
will become less attractive because it will turn out that 
inflation alone cannot describe the observable structure 
of the universe and its dynamics. 

The difficulties that the inflation cosmology faced in 
explaining the current accelerated expansion of the uni- 
verse yielded new cosmological scenarios, such as the 
ekpyrotic [6] and the cyclic [7]. These scenarios arose 
from a set of ideas about "brane worlds" [8], which in 
turn were initiated by research in the still hypothetical 
M-theory. Unfortunately, research shows that brane cos- 
mologies encounter numerous problems, the best solu- 
tion for which is the hypothesis of inflation on the visible 
brane [9]. Thus, these models turn out to be no more 
than exotic variants of the inflation cosmology. 

It is clear from the above that inflation is still the 
best (i.e., the most natural and the most economical) 
hypothesis that explains most, if not all, peculiarities of 
the observable universe. It is therefore sensible to pre- 
serve this hypothesis without introducing additional as- 
sumptions about the existence of a quintessence, which 
would greatly weaken the status of inflation cosmology 
as a descriptive physical theory. U For this, it is reason- 
able to introduce the hypothetical inflaton as the cause 
of the observed accelerated expansion of the universe. 
We note that this idea is unpopular among cosmolo- 
gists. Nonetheless, from our standpoint, the proof of 
the nonequivalence of the inflaton and the " dark energy" 
greatly undermines confidence in the whole inflation sce- 
nario. Indeed, there are no experimental justifications for 
the existence of the scalar inflaton, except for the very 
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existence of a homogenous, isotropic, and practically flat 
universe. Until recently, it was thought that inflation 
solves all major problems of cosmology, and this lent cre- 
dence to the inflation scenario. If it turns out that certain 
global properties of the universe cannot be explained in 
terms of inflation (for example, the current accelerated 
expansion), then we greatly weaken the only (and indi- 
rect) argument that supports inflation (see [10] for an 
excellent survey of the conceptual problems of cosmol- 
ogy). 



II. EXACTLY SOLVABLE COSMOLOGICAL 
MODELS 

An excellent approach to the explanation of acceler- 
ated expansion in terms of inflation was given in [11], 
[12]. The offered models admit a nonpositive deflnite 
self-action potential. This has an interesting implication: 
even if the density of matter in the universe is exactly 
equal to the critical density, a collapse phase may follow 
an expansion phase, whereas expansion lasts infinitely 
long if the potential is positive deflnite. In the presence 
of a negative minimum, the current phase of accelerated 
expansion occurs under a set of natural assumptions [11], 
[12]. 

The developed extended gauge theories of supergravity 
with N =2, 4, 8 allow hoping that the described models 
can have physical meaning. The de Sitter solutions corre- 
spond to the extreme points of the effective potential, and 
the squares of the masses of the scalar fields m for N > 2 
are quantized in units of the Hubble constant Hq. If 
m^/ifg = 0(1), then it is possible to describe the current 
accelerated expansion of the universe [12]. These scalar 
fields have exceedingly small mass (m ~ 10~^^GeV) and 
can significantly change the value of the cosmological con- 
stant during the dustlike ("dark") matter era. Simulta- 
neously, they do not alter the standard predictions of the 
inflation theory, because these fields are far from the po- 
tential minimum in the early universe and " turn on" only 
when the Hubble parameter becomes of the order of |m|. 
Finally, calculations show that quantum corrections to m 
are extremely small. 

Thus, as we have seen, inflation cosmology is rather 
viable. We note that a general study of the dynamics of 
the universe for a given self-action potential is an excep- 
tionally difficult mathematical problem. Therefore, much 
research into model cosmologies where the equations can 
be solved exactly was done during the last seven years 
[13]. The corresponding theories are called exactly solv- 
able cosmological models (ESCMs) or simply exact cos- 
mologies [13]. All ESCMs are based on the extraordinar- 
ily wide gauge arbitrariness of the Fricdmann-Lemaitre- 
Robertson- Walker (FLRW) metric. If we consider the 
general case of universes with such a metric and a self- 
acting scalar field, then by setting the evolution of the 
scale parameter, for example, we can calculate the self- 
action potential that would lead to such an evolution. In 



other cases, the dynamics of the field or of the Hubble pa- 
rameter were fixed. In the latter case, it was convenient 
to introduce a new independent variable, the number e of 
expansions of the universe. We note that the apparently 
boundless arbitrariness is narrowed by the energy condi- 
tions: the weak, the strong, or the dominant condition, 
depending on the context. Therefore, not every possible 
evolution of the scale parameter, for instance, can occur. 
We study this question in greater detail. 

An obvious drawback of ESCMs are "pathological" 
self-action potentials that most likely cannot be justi- 
fied by elementary particle physics. Nonetheless, the 
study of ESCMs can be very important for inflation cos- 
mology. Indeed, the chaotic inflation theory assumes 
that practically any potential satisfying certain condi- 
tions (see the introduction) leads to inflation with exit, 
oscillations, transition to the Friedmann phase of radia- 
tion dominance, and so on. If this is indeed so, then the 
pathological form of the potentials normally appearing 
in ESCMs is insigniflcant. We can gather statistics over 
many ESCMs and make a well-founded claim about the 
reliability of the claims of the chaotic inflation theory. 

The investigations conducted in the works cited above 
led to the following conclusions: 

1. The idea of slow rolling is true. Inflation does indeed 
occur under an extremely broad range of self-acting po- 
tentials, and there is hence no need to fix a certain form 
of the potential to obtain an inflationary universe. ^] 

2. The exit from inflation turned out to be a difficult prob- 
lem. For many model potentials, the universe never stops 
infiating. The exit is generally achieved by fine-tuning or, 
to put it simply, by parameter fitting. 

3. The origin of oscillations is unclear. Even though exact 
models of infiation have been constructed, they definitely 
do not contain an oscillation phase (unlike nonintegrable 
models such as ^ = m^0^/2, where the presence of oscil- 
lations is confirmed by qualitative estimates and numer- 
ical integration [3], [14]). This does not mean that such 
a state cannot be described by any exact cosmology. For 
example, a potential that describes damped oscillations 
of a scalar field was constructed in [15]. But it is doubt- 
ful that such a potential can be obtained from quantum 
field theory or string theory. 

4. Every ESCM with exit has the following drawback: 
after the exit from infiation, the universe cither en- 
ters a radiation-dominated stage, where it stays forever, 
or enters a dustlike matter stage that is succeeded by 
a radiation-dominated stage, whereas it should be the 
other way around. This is a fairly common phenomenon. 
For example, an ESCM with a complex scalar field was 
studied in [16], and the same problem occured.The most 
successful ESCM model to date, described in [15], in- 
troduced an additional hypothesis of the existence of a 
certain phase transition that would deliver the universe 
from a radiation-dominated era to an era dominated by 
dustlike matter. But a physical theory can hardly be 
based on such ad hoc hypotheses. The last problem can 
seem somewhat artificial. Indeed, according to the classi- 
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cal paradigm, the transition from the radiation era to the 
dusthke matter era is normally described by the Fried- 
mann equations, the main contributions to the energy- 
momentum tensor being from the electromagnetic field 
(during the radiation era) and later, when the tempera- 
ture has fallen enough, from nonrelativistic matter (hy- 
drogen). But the difhculties that arise in item 4 have 
to do with the solutions of the equations describing the 
scalar field and gravity, and the scalar field during the 
decoupling stage is very small compared to the radiation 
and baryons. Therefore, its dynamics should not play a 
significant role. 

Modern observational results do not agree with this 
conclusion. It is now known that the main contribution 
to the energy-momentum tensor in the current era comes 
from "dark energy." The observed microlensing shows 
that regular baryonic matter is the next in distribution 
after dark matter. In other words, even at the current 
stage, the dynamics of the universe is governed not by 
baryons but by "dark energy." 

What is this " dark energy" ? A reasonable assumption 
is that it is described by a scalar field of unknown na- 
ture. It is also reasonable to assume that "dark energy" 
has always played a dominant role in the evolution of the 
universe. In other words, not just inflation but indeed the 
entire post inflation dynamics is governed a scalar field 
(or scalar fields). Taking our limited knowledge of the 
universe into account, we naturally attempt to describe 
cosmological evolution using only one scalar field. Be- 
low, we discuss this specific problem, Specifically, we 
propose a simple method for constructing self-acting po- 
tentials that lead to inflation (or to several inflations) . In 
this approach, inflation with exit versions arise without 
parameter fitting of any kind, containing an asymptotic 
Freidmann era characteristic of dustlike matter. 

This method has already been used in the theory of 
nonlinear integrable equations and is known as the Dar- 
boux transformation (DT) [19]-[21]. To avoid misun- 
derstandings, we stress that the goal of this paper is not 
to construct an adequate cosmological model or to pro- 
pose a hypothesis about the physical nature of " dark en- 
ergy" (as is done, e.g., in [12]) but to describe a mathe- 
matical method, which, as we hope, can be used to con- 
struct such a model in the future, if this is indeed possi- 
ble. 



III. FORMULATION OF THE PROBLEM 

We consider the simplest case of a flat universe with 
the FLRW metric 

ds^ = dt^ - a{tf {dx^ + dy^ + dz^) , 

filled with a scalar field with the Lagrange function 



Hence forth, we use the system of units SttG = c = 1, 
commonly accepted in cosmology. We limit ourselves to 
the flat case (k = 0) because, on one hand, this is the 
simplest case and, on the other hand, the curvature in 
the early universe is negligible. It is easy to see that if 
(j) — (j){t), then Einstein's equations with a cosmological 
term 

reduce to the Schrodinger equation 

'ip = 3{V + A) (1) 

where tfj = a^. We assume here that V — V{cj){t)) = V{t) 
and that the time dependence of the field is given by 

^--ll^log* (2) 

In this language, it is very easy to formulate a naive 
method for constructing ESCMs. For example, we assign 
iplt); substituting in |(2J), we determine the time depen- 
dence of the field (f>{t)-, and using (Q) and l|2Jl together, 
we find as a function of (/) in a parametric form. 

The given method looks very simple. But it should not 
be thought that such a potential can be found for any 
evolution of the scale parameter (or the Hubble parame- 
ter H = d/a). First, as is clear from a naive choice 
of a{t) (i.e., ip) most likely leads to situations where cj) 
becomes imaginary after some time. Second, only in ex- 
ceptional cases is it possible to reconstruct V as an ex- 
plicit function of 0, which of course hampers attempts to 
find a physical meaning of such a potential. But in this 
paper, we concentrate on the first point and disregard 
the second. As we see below, choosing (p such that the 
right-hand side of idJ changes sign in a given interval, we 
can construct cosmological models with inflation and an 
inevitable exit from it. 

In what follows, we require the energy conditions. 
Following [22], for example, we define them as 

1. the weak energy condition, p > 0, p + p > 0] 

2. the strong energy condition, p + p > 0, 3p + p > 0; 

3. the dominant energy condition, p > 0, p > p > —p, 
where p and p are the respective density and pressure. 
It is obvious that the weak energy condition (the most 
important one) holds if and only if the right-hand side of 
^ is positive. This simple property is very importance 
later. 

We apply the DT [20], [21] to We have two justifi- 
cations for this. First, the DT is an isospectral symmetry 
of the Schrodinger equation because the DT preserves 
the asymptotic characteristics of the transformed func- 
tions up to a constant multiple. In particular, square- 
integrable functions stay in that space if the support 
functions that generate the DT do not have zeroes or 
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FIG. 1: 



have alternating zeroes if the number of functions is even 
[23]. This property can be used as follows. Let the 
scale parameter have an asymptotic characteristic of a 
flat universe filled with dustlike matter, i.e., a ~ t^/^. 
Dressing the appropriate solution of Eq.(^ with a sup- 
port function from L2 (any regular solution of with 
the same potential suffices, even if it does not belong 
to L2), we obtain a new solution of Einstein's equations 
with the same asymptotic condition. In other words, 
even though the behavior of the scale parameter can be 
rather odd (for example, it can grow rapidly on a finite 
interval, and the universe therefore experiences an infla- 
tionary period), the universe enters a Friedmann expan- 
sion era with the two-thirds law in any case. footnoteOf 
course, the support function must have no zeros between 
the inflationary era, where a > ,and the Friedmann era 
proper. Otherwise, the inflationary era and the matter- 
dominated era are separated by a singularity, i.e., repre- 
sent two casually unconnected universes. This property 
guarantees an exit from inflation. 

Second, DTs can be used to construct inflation models. 
For this, we rewrite Eq.(|2Jl in the form (fp' = F^"'^' — V, 
where V^^^ is the potential constructed from the poten- 
tial V using Crum's formulas and the support function 

(see formulas © below) [21]. We assume that the 
function 7p has N zeroes. Then V^^'> has at least TV sin- 
gularities. We assume that these singularities are not 
present in the initial potential V . If V'^^'^ is not positive 
definite in the neighborhood of one of these singularities 
(or V has a singularity of its own, around which it is 
positive definite), then there exist ti and t2 on different 
sides of the singularity such that ^(ti) = </>(i2) = 0, for 
the (^^(i) < 0, ti < t < t2. The proof of this statement 
is obvious. It is now clear that if the given conditions are 
met, then the universe certainly enters an inflationary 
phase when t > ti. 

Indeed, we assume that the universe is born at t — ti 
{U < ti) and that there exists a period of time during 



which all energy conditions are met. It can be shown 
that as we approach ti, the strong energy condition is 
violated first, while the weak one holds. As a result, the 
universe expands up to the instant ti. At t — ti, we 
obtain the de Sitter phase, but for only one instant. A 
forbidden zone then follows, which should be removed 
by matching the solutions at t = ti and t = t2- We note 
that this matching requires stricter conditions than those 
found in corresponding problems in nonrelativistic quan- 
tum mechanics based on the Schrodinger equation. Here, 
we require that not only the wave functions and their first 
derivatives but also the second derivatives be the same 
at the matching point. The latter is necessary to avoid 
discontinuities of density and pressure. The general so- 
lution of (^3) has two constants of integration; therefore, 
it is not obvious that we can always satisfy all three con- 
ditions. But it is easy to verify that the third condition 
(equality of the second derivatives) follows automatically 
from the first two! This fortunate circumstance is a con- 
sequence of matching the solution at the inflection point 
of log'f/'. The price we pay (a discontinuity of the third 
derivative of the scale parameter) does not seem exces- 
sive, at least while we consider the system as a classical 
(i.e., not quantum) system. 

Thus, controlling the zeroes of the support function 
of the DT, we can construct exact cosmologies with in- 
fiation. If solutions with known asymptotic characteris- 
tics (for example, evolving according to the two-thirds 
law) are used as initial base solutions, then the universe 
necessarily exits inflation, without any parameter fitting. 
Below, we demonstrate examples of how this approach 
works. 



IV. A SIMPLE MODEL 

The law a ~ t^^^ gives the potential 



The general solution of 1^ with potential ^ and without 
a cosmological constant has the form 

^(-) = -^-c2^^ (4) 

where ci and C2 are the constants of integration. The 
solution describes two disconnected universes (Fig.l). 
The first is infinite in the past but has a singularity at 
to — — (ci/c2)^/'^. The second is rather strange. It is 
born from a singularity at i = tg, then expands, and un- 
dergoes superinfiation, i.e., the scale parameter becomes 
infinite at t = 0. Then the universe contracts, the scale 
parameter changing sign. Of course, this solution cannot 
be taken literally. The superinfiation zone (to the left 
and right of t = 0) is forbidden because the weak energy 
condition is violated in it. It is easy to flnd the density 
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and pressure 

Ci (ci + 4c2t^) 



P 



(Cl + C2t3f ' ■ 3t2 (ci + 02^3)' 

Therefore, the above mentioned forbidden zone is in the 
interval between ti and ^2, where 

T 1/3 / \ 1/3 



P ■ 



Cl 



V2 



-0.5 



^2 



2+4 

V2 



1/3' 



1.6 

C2 



1/3 



(5) 



We consider solution Q to hold only for t < ii, where ti is 
given by Q . To continue the solution into the region t > 
<i, we rewrite Q with different constants of integration, 
denoted by 6i and 62: 



h 
t 



(6) 



For (6), the weak energy condition is violated in the inter- 
val (^'1,^2) where the t'^ 2 ^'^^ determined by expressions 
with the substitutions ci — > 61 and C2 &2- The 
constants hi and 62 should be chosen such that they sat- 
isfy H 

i,^-\t{)^i^(+\t'2), i^(-\h)=i^(+\t'2). (7) 
It is noteworthy that another relation, 

4,(-\h)^i^(+\Q. (8) 

is a consequence of ©. It thus follows that when match- 
ing ip'^~\ti) and (^2), the density, pressure, and rate 
of density variation are continuous! 

An obvious deficiency of this solution is a jump in time: 
ti ^ t^- Therefore, the solution is not defined on the 
interval {tx^t'^)- But it is easy to move t in and 
such that the two conditions that 

1. the initial singularity corresponds to time t = and 

2. t\ — t'2 with no jump in time. 

are met. We skip the calculations and present only the 
final answer: 

1. For < t < ti, the dynamics are described by the 
formulas 



Cl 



t - Ci/C2 



C2 {t - C1/C2) 



2/3 



(9) 



{t-c,/c2r 



2. For ti < t, the dynamics are described by the formulas 




Here, 
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Cl 



= — 1 + 2-f ^ 



C2 



V2 



1/3N 



(? = — f 1 + f 16 + 12V2^ 

C2 



2.6^, 

C2 



4.2^. 



(11) 



C2 



we introduce the new constants cf 2 for convenience, and 
t2 = ti .The matching here is 



^Mit) ^ ^ ^3(3 + 2V2) (t - g)\ p'^^Hh) - P^+Hh), P^-Hh) = P^^Kh). 

(10) The graph of a as a function of t for the second universe 

T/(+) _ is given in Fig. 2. The curve is described by the formula 

~ (i - gf ' a(-) = V'^'HO for t < ti and by a(+) = ^^+\t) for t > 
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ii, where i^'-~\t) and V^+H*) 

are given by ij^ll and (|10|l . 

We specifically choose ci = C2 — 1. Figure 3 describes 
the acceleration a as a function of time. The break is at 
t = ti, given by l(TT|) . 

As can be seen, inflation occurs on a finite interval 
(where there is a tooth on the t axis) .We also note that 
as t cx), the potential behaves as 
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We have thus constructed a simple model of inflation 
with exit and an asymptotic transition to the Friedmann 
expansion era, located in the matter-dominated phase. 
Of course, this example does not describe our universe, 
because there is no oscillation phase. Neither is there sec- 
ondary acceleration. But it is evident that this method 
can be used without fundamental difficulties to construct 
model cosmologies with exit without parameter fitting. 
Neither it is necessary to use hybrid models. We note 
that the solution describes another, collapsing universe 
with no beginning but with an end. These solutions 
seem rather strange, but their appearance is inevitable 
because the equations of general relativity are invariant 
under time reversal. 



V. THE DARBOUX TRANSFORMATION 

It is clear that any exact solution of (QJ satisfying the 
energy conditions generates an exact solution of the Ein- 
stein equations in the FLRW metric. Because not all 
solutions i/j of generate a positive definite expression 
after substitution in the right-hand side of the exact 
solutions of the Einstein equations are a subset of the ex- 
act solutions of . On the other hand, practically all ex- 
act solutions of Schrodinger equation with " physical" 
potentials can be constructed using a DT J^. Examples 
are the harmonic oscillator, multisoliton, and even finite- 
gap potentials. We use this as an "experimental fact" 
without attempting to give precise mathematical mean- 
ing to statements like "all potentials for which Eq.QJ 
can be explicitly solved for all values of the spectral pa- 
rameter are constructed using the DT and, possibly, an 
additional finite set of symmetries from other potentials 
with this property." Informally speaking, we can then 
conclude that all exact solutions of the Einstein equa- 
tions in a flat FRLW universe filled with a real scalar 
field can be constructed using DTs. 

As a simple example, we consider potential which 
is constructed using a single DT from a zero potential 
1^(0) = 0. The result of a DT applied N times can be 
written compactly using Crum's formulas: 



where 
A^v = det 



A 



Af+l 



A 



N) 



V - 



N 



2 In An 

3 ~ 



k = 0,...,7V- 1; 



dt^ 



z = l,...,7V, 



and 



Ajv-i-i — det 



dSk 
df" 



k = 0,...,7V; 



l,...,iV+l. 



In these expressions, the ipi are particular solutions of 
with a common initial potential V and different eigenval- 
ues 3 A. It is now clear how we can use DTs to construct 
models with two accelerations. We already mentioned 
that the potential y'^^^must have singularities. It fol- 
lows from Crum's formulas that all singularities coincide 
with the zeroes of functions that are contained in Crum's 
determinant An Therefore, the problem is reduced to 
controlling these zeroes. We assume that the initial po- 
tential is given by lj2Jl. Solution 10} corresponds to a 
trivial cosmological term. The solutions of Q with pos- 
itive (A — k"^ /i) and negative (A = — A;^/3) values of the 
lambda term are 



^4 



1 



1 

ki 



kt 



M2 1 



1 

ki 



-kt 



sin kt 
kt 



cos kt 



V2 



cos kt 
kt 



sin kt 



(12) 

where /ii_2 and 1/1^2 are arbitrary constants. Using these 
functions for different values of k in formulas (|V| and 
regarding -iAN+i defined by Q for all N , we can con- 
struct a set of solutions with the Friedmann two-thirds 
asymptotic behavior for t —t +00. 

It is clear that the functions ip'^'^^ contain any num- 
ber of singularities near which the weak energy condition 
is violated and acceleration therefore occurs. As in the 
example in Sec. 3, the solutions must be cut off at the in- 
flection points Mj of logV'*-^\ where the exact de Sitter 
phase is realized, and matched along these points. We do 
not need to match solutions at adjacent points. By adja- 
cent points, we mean points Mj and Mj^i between which 
there is only one singularity of As shown in Sec. 3, 

we can always match any two points Mj , which gives a 
degree freedom in constructing the required solution. 

We note that the explicit form of the solution quickly 
becomes very complicated as N increases. We consider 
some examples. Henceforth, we suppose that ci = C2 = 1 
everywhere. 

A single DT with a support function tp+ from l(T^ 

gives 



{x-iy + k^ + i 



{x + ly + k^ -1 



2fc2 ^2g2 



^{x — l)e^ + (x + l)e 
- e"^^ - 4^x^ - 2^ 



3 {^{x ~ l)e'' + {x + l)e-^) 



2 ' 



where ^jli = 1, 112 — IJ^, x — kt. For /x < is given in 
Fig. 4. Here, we already removed the regions where the 
weak energy condition is violated. The boundary points 
(where (/> = 0) must be matched. It is then clear that 
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such a solution describes two universes. The first con- 
tracts to a final singularity, experiencing a time-reversed 
inflation at a certain instant prior to the final collapse. 
The dynamics of the second universe is a time reversal 
of the evolution of the first: it is born from a singularity, 
expands for some time, then undergoes inflation, exits, 
and settles asymptotically on the Friedmann expansion, 
characteristic of a matter-dominated universe. 

A new situation arises for fi > 0. Figure 5 shows the 
graph of V'^^H^) for /u = 2. It demonstrates the evolution 
of three universes. The first collapses into a finite singu- 
larity but has no beginning. The second is born from an 
initial singularity, expands for some time, then begins to 
contract, and collapses into a final singularity (we recall 
that k = 0). The third is born, expands, undergoes infla- 
tion (the boundary points are matched), and exits into a 
two-thirds phase. 

As a final example, we present a solution obtained by a 
double DT (Fig. 6). One support function is a solution of 
(1) with A = +1/3, /ii = 1, /i2 = 2 and k = 1; the other 
has A = +4/3, fii — 1, fj,2 — 3 and k — 2. We also have 
three universes here. We note that the second exists for a 
finite period of time, during which it undergoes inflation 
and then collapses. These solutions are very similar to 
the solutions in [11] that we mentioned at the beginning 
of Sec.2. 

We have not attempted to construct models with two 
(or more) accelerations. It is clear that this problem can 
be solved in principle using this method. Moreover, the 
acceleration rate is determined by the density p; there- 
fore, constructing potentials with several singularities, we 
can obtain several inflations, each with its own intensity. 
If the second inflation corresponds to a smaller value of 
p than the first, then it is much slower. But we think 
that it is too early to begin constructing models using 
DTs. It is necessary to understand how an intermediate 
oscillation phase can be included. It is still unclear to 
us how can this be done. It is obvious that any inflation 
model without such a phase is impracticable. |9j 
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VI. CONCLUSION 



We suppose that the exact solution a{t) — ^•\/-0(^)(i) 
describing two acceleration phases and an intermediate 
oscillation phase is constructed. We obviously encounter 
the problem of interpreting the potential y (-'^^ , which is 
undoubtedly hideous. It is unlikely that it can ever be 
described by any field theory or even string theory. In- 
deed, even the potentials obtained by a single DT are 
already too complicated and unphysical, and they are 
simple enough to describe even two consecutive acceler- 
ations without oscillations! 

A natural reaction to the appearance of such expres- 
sions is the conclusion that real-life cosmology is not ex- 
actly solvable. From this standpoint, exact cosmologies 
can be useful models, but they are incapable of describ- 
ing the real universe. But just how correct is this con- 
clusion? It is quite possible that the correct potential 
of a self-acting scalar field is not defined by elementary 
particle physics but appears for a different reason. Any 
suppositions of this kind are speculative and unreliable. 
Nonetheless, we attempt to point out a source of such 
strange potentials outside both field theory and string 
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theory (of the future). 

According to the chaotic inflation theory of Linde, all 
possible vacuum states are realized in different areas of 
the universe, which in turn must lead to different mass 
spectra of the elementary particles populating these re- 
gions. In this paradigm, our region of the universe cannot 
be identified using physical laws alone; a certain "weak 
anthropic principle" is also required. In other words, even 
if an all-encompassing M -theory can be constructed and 
it might justify our expectations, it would still prove in- 
sufficient to explain why the vacuum state of our universe 
is the way it is and not different, and so on. 

To illustrate this picture, Linde introduces the con- 
cept of the multiuniverse[14]. In the simplest case, the 
multiuniverse can be described as an infinite sum over 
all possible actions describing all possible quantum field 
theories and M -theories. ^3 Even though hypotheses of 
this kind can seem speculative, it must be borne in mind 
that all existing explanations of the observed density of 
the vacuum energy 0.7pc) have been obtained using 
the anthropic principle, which apparently makes sense 
only if the multiuniverse hypothesis is true. [TH 

If we accept this hypothesis, then it is clear that the 
existence of regions with the strangest self-acting poten- 
tials is possible. In this case, we must readdress the ba- 
sic questions we normally ask. For example, potentials 
V{(j)) are generally taken from particle theory {V ~ 0^ 
and V ^ (j)'^) or string theory V ^ e""^. Why these and 
not V ^ (j)^, for example? Because a theory with such 
a potential is nonrenormalizable for d = 4. We have ev- 
ery reason to believe that a nonrenormalizable theory is 
a "bad" theory because it requires a new Lagrangian on 
each new scale and therefore contains an infinite num- 
ber of indefinite parameters as a rule. As well as being 
renormalizable, models must also satisfy other generally 
known properties: unitary, a Hamiltonian bounded from 
below, etc. 

We now suppose that we have succeeded in construct- 
ing a family of potentials giving an evolution of the uni- 
verse consistent with observations (for example, using 
DTs). The next step would be to check these models 
using the methods of quantum theory (renormalizability, 
etc.). If one of these potentials satisfies the appropriate 
checks, then we can perfectly well use it even though we 
did not obtain it using the basic principles of field theory 
or string theory. We obtained it using the multiuniverse 
concept, and field theory (or string theory) is only re- 
quired to verify its consistency. 
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[1] The alternatives to a quintessence are models with a 
nonzero cosmological term. Unfortunately, in terms of 
modern field theory(or string theory), it is impossible to 
obtain a model with such a small vacuum energy density 
that would not confiict with observations 

[2] This was done in Guth's early models, where the inflation 
was tied to the Higgs form of the potential[l]. 

[3] From this standpoint, such elegant models as Chaply- 
gin's gas model[17] are unsatisfactory, if not necessarily 
untrue. The same can be said of hybrid models[18]. 

[4] In the context of nonlinear integrable equations, it is 
proper to call it the Darboux-Crum-Matveev method. 

[5] The existence of a "forbidden" zone correlates with the 
known lack of superinfiation in cosmologies with scalar 
fields 

[6] It is easy to verify that O does not hold if bi = ci and 
62 = C2; therefore, we must introduce the second solution, 



Eq.linj. 

[7] We do not consider potentials such as the rectangular 
well, limiting ourselves to continuous functions. 

[8] A common idea for constructing many exactly solvable 
cases is form-invariance. For example, the harmonic oscil- 
lator and finite-gap potentials are constructed by closing 
a chain of discrete symmetries. The same is true for cer- 
tain Painleve transcendents. 

[9] A certain hint is given by the astonishing potentials ob- 
tained by closing the Veselov- Shabat chain[24]. These 
potentials contain amplitude-modulated oscillations and 
must be considered exactly solvable. The spectra com- 
prise a set of several arithmetic progressions, and de- 
pending on a special condition, the solutions either can 
be expressed in terms of Riemann's function or are gen- 
eralizations of the Painleve transcendents. In particular, 
the solution of a three-hnk (JV = 3) chain can be ex- 
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pressed in terms of the transcendents of the Piv equa- 
tion. 

[10] In fact, Lindo goes even further, proposing to include 
even those models with no description in the Lagrangian 
formalism and even discrete models such as cellular au- 
tomata. Some do not stop even here, considering all pos- 



sible mathematical structures to be elements of the mul- 
tiuniverse(see[26]). 
[11] An unusual alternative is the "self-consistency" hypoth- 
esis proposed in [27]. 



